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^ ! Abstract: A system of spinless fermions in d = 1 + e dimensions, at zero-temperature 

^1 and in random external potential is studied using the perturbative renormalization group 

c/3 I to first order in disorder and to second order in interaction. We find a superconductor-to- 

Anderson insulator quantum fixed point at an infinitesimal value of disorder and calculate the 
ri . correlation length and the dynamical exponents to the lowest order in e and in interaction. 

The scaling of conductivity with temperature and the behavior of characteristic temperature 
scales on both sides of the transition is determined. The model may have relevance for a 



O ' p-wave superconductor at low temperatures in strongly disordered media. 
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The effect of Anderson localization on superconductivity or superfluidity represents a 
fundamental but still incompletely understood problem in condensed matter physics. The 
effects of disorder may be expected to be particularly pronounced in low dimensions, d <2, 
where arbitrarily weak disorder localizes all single-particle states |l[] in a non-interacting 
system. Experiments on thin, microscopically disordered films at low temperatures, for ex- 
ample, have revealed a sharp separation between the insulating and superconducting regimes 
as some parameter controlling the amount of disorder is varied [§ . These measurements are 
usually interpreted as consequences of the proximity of a quantum (T = 0), superconductor- 
insulator (SI) critical point f^, [§]. Properties of such a critical point have been a subject 
of intense investigations in recent years 0, particularly in the context of the theory of 
superfiuid-Bose glass transition, which assumes that the important degrees of freedom close 
to the transition are bosonic, preformed Cooper pairs. One of the obstacles in better quan- 
titative understanding of the quantum SI transition lies in the fact that the assumed SI 
fixed point resides in the region of strong disorder, and is not easily revealed in perturba- 
tive calculations [Q, 0. A notable exemption from this rule is the SI transition in one 
dimension [0) I0- At T = and without disorder (i = 1 is the lower critical dimension 
for superconductivity; only power-law superconducting correlations are possible since the 
true long-range superconducting order is eliminated according to Coleman-Mermin- Wagner 
theorem. In the language of renormalization group, this is manifested in scaling towards 
finite attractive fixed points . When disorder is added, an immediate consequence of the 
special situation in d = 1 is that the SI fixed point is located at zero-disorder (and at finite 
interaction), suggesting an approach to the problem of the SI transition in higher dimensions 
via expansion around d = 1. Such an expansion would be of particular value since there may 
not be a finite upper critical dimension for the problem P] . 

In this Letter we address some of the above issues by studying a model of interacting, 
disordered spinless fermions in rf = 1 -|- e. The model may be appropriate for the description 
of the SI transition from a superconductor state with p-wave symmetry. Without disorder 
and for e > 0, the interaction coupling constant is irrelevant in infrared if positive, and runs 
towards — cx) if negative, indicating the formation of superconducting long-range order. In 
the presence of disorder, we find a fixed point that controls the quantum SI transition at 
the value of dimensionless resistance of order e, and of attractive interaction of order unity 
(see Figure 1). The critical exponents that characterize the fixed point are evaluated to 
the lowest order in e and perturbatively in strength of interaction. The correlation length 
exponent diverges as e — >■ 0, and the dynamical exponent acquires a positive correction of 
order e. The fact that the fixed point is infinitesimal in disorder enables one to obtain the 
low-temperature scaling of conductivity with temperature in simple way. Due to the pair- 
breaking nature of the disorder in the problem, superconducting gap vanishes at the critical 
point. 

We consider a system of spinless fermions in (i-dimensions, 1 < c? < 2, in random 
external potential distributed according to the probability: 

P[V{f)] oc exp - / d'^rd'^r'V{r)W-\f- r')V{r'). (1) 



The quantity of interest is the zero-temperature free energy G = —{InZ), where brackets 



denote averaging over disorder configurations. Utihzing the standard rephca trick [T^ , after 
averaging over disorder we arrive at the zero-temperature action for fermions: 

S = Sq + Dint + Odis (2) 

where, 



1 J Zn 



^ - ^ d'^Kduj 



and 

TV 4 rl rl 

a,/3=l"' m=l l^^J 

Frequency integrals run over the whole real axis and the measure in i^-space in d-dimensions 
is defined as: 

d'K=—jzTr dk {sin ey-'de, (3) 

where e{K) = k = \K\ — Kp, T{x) is the factorial function, and we assume a circular Fermi 
surface and the ultraviolet cutoff A << Kp. We set Kp = m = h = 1. Greek indices 
enumerate replicas and the zero-temperature free energy is G = /^(exp(— S") — l)/N, in the 
limit of vanishing number of replicas N ^ 0. Note that the individual frequencies, and 
not only their sum, are conserved in Sdis, since the scattering over random impurities is 
elastic. This results in coupling function W having an engineering scaling like frequency (or 
momentum) while F is dimensionless ||TT| , which indicates a finite scattering rate of the plane 
waves. We assume that the incoming (and consequently, the outgoing) momenta in both Sdis 
and Sint are always opposite in direction, since precisely these interactions are responsible 
for superconductivity [|ll| and weak- localization |l[], |jl2[, when present separately. Both F 
and W essentially depend only on angle 6 between the momenta in their arguments, due to 
the restriction on all momenta to the vicinity of the Fermi surface, and F{6) = —F{9 + vr). 
When (i — i> 1, the theory becomes one-dimensional since / d'^K — * j_p^{dk /2tx) J d9{S{9) + 
5{6 — vr)). At (i = 1, the interaction is described by a single coupling constant g = 2{F{6 = 
0) — F{9 = 7r))/7r, while disorder is determined by two coupling constants: Wb = 2W{9 = n) 
for the backward, and Wf = 2W{6 = 0) for the forward scattering over the random potential. 



Hereafter we set Wf = 0, since forward scattering is unrelated to localization in d = 1 |]T3| 
and drop the subscript "b" in Wb- In rf = 1 + e, < e << 1, for F{9) and W{9) smooth 
around 6 = ir the theory is still parametrized with only two coupling constants g and w, 
as defined above. To the lowest order in e and disorder, /3-functions then assume a general 
form: 

/?, ^ ^ = -eA{g) + wB{g) + 0{e', w\ ew), (4) 

P^^^ = wC{g) + 0{e^w^ew), (5) 

d^ 

where w = [w/k) > is the dimensionless disorder- variable, and ^ = ln(A//c) = ln(A/ci;), 
where k (uj) measures the momentum (frequency) of the external legs in the renormalized 
disorder (interaction) vertex. We calculate the requisite functions using the standard dia- 
grammatic perturbation theory in c? = 1, to the second order in g: 

A{g) = '-^ + 0{g% (6) 

B{g) = -^g + 0.10g' + O{g'), (7) 

C{g) = l + g-^g' + Oig'). (8) 

Function A{g) is positive, so at u) = 0, in the infrared limit [C, —>■ oo) g -^ — cxd if initially 
negative, indicating the opening of a gap in spectrum and the long-range superconducting 
order at T = 0. The term proportional to e in (3g comes from the increase of the phase 
space available for the particle-particle, relative to the particle-hole, scattering when e > 0, 
and not, as usually, from g acquiring a finite engineering dimension. Recall that in two 
and three dimensions diagrams of particle-hole type are suppressed by powers of A/Kp due 
to the restriction that all scattering processes happen in the vicinity of the Fermi surface 
|TT|] . In d = 1 the Fermi surface is reduced to two points, and one- loop particle-particle and 



particle-hole diagrams are precisely equal and opposite in sign, so they cancel out exactly. 
This cancelation persists to all orders in perturbation theory [Q, so /?g = in c? = 1 
without disorder. For infinitesimal e however, the homogenous system at T = may be 
expected to behaves qualitatively like in higher dimensions, and our lowest-order result for 
A{g) conforms to this expectation. The second term in f3g is positive for ^f < 0, up to 
two-loops for B{g). This is a pair-breaking property of the disorder, which therefore always 
inhibits superconductivity in this model. Also, C{g) becomes reduced for g < 0, indicating 
a suppression of localization by the attractive interaction in turn. The disorder-variable w 
scales precisely like the dimensionless resistance in li = 1, and the eq. 5 for g = reduces 
to the statement of the Ohms law. In d = 1 + e the resistance and w should be expected 
to be related as i? oc wk'^, but this distinction between the two variables is of order of e^ in 
the /3-functions. Staying infinitesimally close to c? = 1 we thus may still identify w as the 
resistance, i. e. the right scaling variable. It is also interesting to note that the second-order 



(two-loop) corrections in B{g) and C{g) produce only small quantitative, but no qualitative 
changes. 

When e = the problem has been studied before in refs. 7 and 8 within the bosonized 
version of the theory (2) via mapping onto the equivalent system of two-dimensional classical 
Coulomb plasma. Our perturbative results are in general agreement with these calculations. 
Some numerical differences originate in different implementation of the renormalization pro- 
cedure: the coefficient of the ^f^-term in C{g) is different, for example. The topology of the 
scaling diagram, the behavior of the conductivity at the transition and the divergence of the 
correlation length are the same as in our analysis at e = 0, these differences notwithstanding. 

The main new feature brought by finite e is the appearance of the finite-disorder fixed 
point in the physical region of the parameter space that determines the critical behavior at 
the quantum SI transition (Figure 1). It is located at 



9c 



1 - \/3, (10) 



where the numbers are obtained on the basis of the perturbative series in eqs. 6-8. Note 
that the fixed point is infinitesimal only in disorder and not in interaction, so our results are 
not exact even to the lowest order in e. Interaction at the fixed-point is of order dr = 1, 
the number of dimensions in directions of imaginary time. A more precise determination 
of the scaling equations would require a calculation and resummation of higher order terms 
in g. Under the assumption that there are no intermediate fixed points, scaling trajectories 
starting above the separatrix on Figure 1 end in the non-interacting, strongly disordered 
sink, which we identify with Anderson insulator. This should be contrasted with the Bose- 
glass phase in which tightly bound pairs of fermions would be localized. For initially positive 
interaction g, the scaling is always towards the insulating phase, unless w = 0, when fiows 
terminate in the trivial, unstable fixed point at the origin. 

The relevant direction at the SI fixed point plays the role analogous to temperature in 
thermal phase transitions. The linearization of the fiow in the vicinity of the fixed point 
yields the correlation length exponent, L^or oc 6~'^: 

(11) 



f^(9c)C'(9o) v^ ' 

where the prime denotes the derivative of the function and the parameter 6 measures the 
distance from the transition. When e — > the exponent u diverges, indicating the vicinity of 
the lower critical dimension for the transition. The value of the exponent satisfies the Harris 
criterion u > 2/d [^ for all values of e. To the lowest order in disorder and interaction, the 
dynamical critical exponent acquires a positive correction: 

z = l + ^ + Oiw,g,,wl), (12) 



as a result of different treatment of frequences and momenta by the disorder term in the 
action (2), which breaks the relativistic invariance of the homogeneous system. 

We may define the superconducting quasi-transition temperature Tgc (proportional to 
the superconducting gap at T = 0), as the frequency (momentum) scale at which g diverges. 
When w = 0, it follows from the eq. 4 that: 

To^sc = A exp — . — — — , (13) 

e|5'o|ln2 

where (70 < is the initial value of g. True superconducting temperature is of course zero 

for d < 2 due to the effect of thermal fluctuations which are not included properly in our 

zero-temperature theory, nevertheless, the above defined quantity sets the temperature scale 

at which the Cooper pairs start to form. On the superconducting side, close to the SI 

transition this energy scales to zero as Tg^ oc 5^^ . The diverging size of the Cooper pairs at 

the transition is a consequence of the pair-breaking nature of the disorder in our model. It 

thus seems unlikely that the SI transition in this systems can be described by an effective 

bosonic theory in terms of preformed Cooper pairs. 

Under the assumption that at low, but finite temperatures, scaling equations do not get 

substantially modified, we may infer the temperature dependence of the conductivity [||. 

At the superconducting side of the transition, w ^ 0, and the conductivity is given by the 

sum of classical (Boltzmann) part and quantum (weak- localization) correction: a = as + crQ- 

According to our calculations w vanishes at some finite temperature Tsc, at which g diverges. 

At this temperature superconducting fluctuations become strong which results in dramatic 

reduction in resistivity. Within our approximation, as actually diverges at Tgc, while ag ~ 

— L^"*^, where Lt oc T~^ is the inelastic scattering length, remains finite. Thus a = as close 

to Tsc, and in our units o"b = l/(wA;^~^), with temperature scaling as momentum T = k. 

The behavior of w close to Tgc is determined by the scaling equations in the region g << —1. 

From the eq. (8) C{g) ~ — 5'^/2 in this region, implying the resistivity, to the lowest order 

in e: 

2T 

P = o-~^ ^ WqT exp — — — 7^^ — ;;r^- (14) 

^ " ^ e2(ln2)2(T-T,J 

For more general A{g) and C{g), we would still expect to find an activated form of the 
resistivity close to Tgc, but with the power of (T — Tgc) and the constants in the exponent 
modified. 

On the insulating side of the transition, at low temperatures the system is in the localized 
phase and the transport proceeds via hopping between the localized states. The conductivity 
then should have an activated form, similarly (but not identically) to the resistivity in the 
superconducting phase close to Tsc'- 

cr~exp-(^)", (15) 

where the power a depends on the hopping mechanism. The temperature scale Tioc ~ l/Lf^^, 
where Lioc is the localization length at the insulating side. This length may be defined as a 



scale at which w becomes of order of unity upon scahng, so Lioc oc Lcor and Tioc also vanishes 
as Tioc (X 5^''. In d = 1, Lioc oc exp(const/5^/^), akin to the divergence of correlation length 
in the Kosterlitz-Thouless transition p. 

Right at the SI transition the resistance w has a finite value Wc oc e. The conductivity 
is still determined by the weak-disorder expression. Now as = T'^^""^^'^ /wc as T — i> 0, and 
the quantum correction aq ~ T^^^"*^)/^ is still negligible infinitesimally close to d = 1. The 
conductivity at the point of SI transition thus diverges as: 

'J^~{1— const. e) 

a = , , (16) 

Wc 

with const. ^ 1.20, to the lowest order in g. When e = 0, we find a = a^ oc (In T)^/T right 
at the transition. 

The results presented in this paper are not directly relevant for the experiments on 
thin films [0 due to the assumed spinless nature of fermions in our model. Disorder is not 
pair-breaking for the usual s-wave pairing, and the SI transition in that case should be in a 
different universality class, for which the theory in terms of the bosonic Cooper pairs becomes 
viable 0. The transition between a p-wave superconductor and insulator, however, could 
be expected to be in the same universality class as the considered model. ^He in strongly 



disordered media |]I6[ such are aerogels may offer an experimental realization of the quantum 
transition we considered. 

In summary, we studied a 1 -|- e-dimensional model of interacting, disordered fermions 
which exhibits a direct, continuous, superconductor-Anderson insulator transition at zero 
temperature. The value of parameter e controls the value of resistance at the quantum fixed 
point, and we calculated the critical exponents to the lowest order in e and interaction. The 
behavior of conductivity at both sides of the transition is obtained, and the superconducting 
gap at the critical point is showed to vanish. The model may have relevance for p-wave 
superconductors in disordered media at low temperatures. 

We acknowledge useful conversations with Professors I. Affleck, P. Stamp and Z. Tesano- 
vic. The author thanks NSERC of Canada and Izaak Walton Killam foundation for financial 
support. 



Figure Captions: 

Figure 1: Schematic scaling diagram in the interaction-disorder plane in d = 1 + e for 
spinless fermions. The scaling trajectory flowing right into the flxed point separates the 
flows towards the superconducting {g = —oo,w = 0) phase from the flows towards the 
non-interacting localized insulator {g = 0,w = oo). 
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